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CYCLIC STRUCTURAL ANALYSES OF AIR-COOLED GAS 


TURBINE BLADES AND VANES 


by Albert Kauitaan and Raymond E. Gaugler 
MSh Levjis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

T!.e creep-fatigue behavior of a fully 
impingement -cooled blade for four cyclic cases 
was analyzed by using the Bias 55, finite- 
element, nonlinear structural computer program. 
Expected cyclic lives were calculated by using 
the method of Strainrange Partitioning for re- 
versed inelastic strains and time fractions for 
ratcheted tensile creep strains. Strainrange 
Partitioning was also applied to previous re- 
sults from a one -dimensional cyclic analysis of 
a film-impingement -cooled vane. The aiialyses 
indicated that Strainrange Partitioning is mor- 
applicable to a constrained airfoil such as the 
film-impingement -cooled vane than to the rel- 
atively unconstrained fully impingement -cooled 
airfoil. 
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IMLASTIC STRAINS can be indueed in airfoils as 
a reeulfc of high themal stresses during engine 
transients and creep during steady-state opera- 
tion. Repetition of inelastic straining under 
normal engine operation results in progressive 
creep-fatigue damage end eventual failure by 
cracking. A major difficulty in applying cyclic 
life prediction methods such as Strainrange 
Partitioning (l to 3)* to a turbine blade is 
that the computation of inelastic strain cycles 
is a formidable \indertaking. 

In advanced aii’Craft engines, the high- 
pressure -stage turbine blades and vanes involve 
geometries and envlroniaents that require the use 
of three-diiaensional, finite -element, stress 
analysis prograitis. Nonlinear, structursQ. anal- 
ysis coirrputer programs are necessarily complex 
and, for large problems, use substantial aitiounts 
of computer tiiae and storage. 

One -dimensional stress analyses based on 
strength-of -materials theory have been, extensiv- 
ely used in the design and cyclic analyses of 
cooled turbine blades because of their relative 
simplicity and because nonlinear three- 
dimensional structwal prograias have not been 
available until recent years. The one- 
dimensional analyses were adequate, provided 
the section analyzed was a chord width away 
from the end and the spanwise metal temperature 
gradients >;ere reasonably linear. In (4) and 
(5), strain hysteresis loops were computed from 
one “dimensional prograsiis for air-cooled vanes 
and blades with airfoil aspect ratios of ap- 
proximately 2, However, a one -dimensional 
stress analysis would be of questionable valid- 
ity for cooled blades having aspect ratios of 
the order of 1 and nonlinear spanwise temper- 
ature gradients, as is discussed in ^6). 

Strainrange Partitioning, which has been 
developed at the MSA Lewis Research Center by 
Manson, Halford, and Hirschberg, is a method 
for calculating low-cycle fatigue life under 
conditions of coiuplete inelastic strain revers- 
al. The inelastic strain hysteresis loops are 
divided into foui* basis components: 

(1) Tensile plastic strain reversed by 
compressive plastic strain, Aepp 

(2) Tensile plastic strain reversed by 
compressive creep strain, ^pc 

(3) Tensile creep strain reversed by com- Kaufman 
pressive plastic strain, Z^e^p 

Numbers in parenthesis designate 2 

References at the end of paper. 
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(4) Tensile creep strain reversed by com- 
pressive creep strain, 

In cases where complete inelastic strain revers- 
al does not take ^ace, additional consideration 
ranst be t^iven to the -unreversed component of 
strain through a ductility exhaustion or tiitse 
fractions approach. 

The purposes of this study were, first, to 
demonstrate the use of an advanced, nonlinear, 
stress analysis computer program in calculating 
blade life and, second, to evaluate the effects 
of airfoil thermal and mechanical loading con- 
ditions and geometry on the inelastic strain 
cycles and the general applicability of the 
Strainrange Partitioning method. It is not the 
intention of this paper to verify the accuracy 
of the Strainrange Partitioning method since 
this is beyond the limited experiinental infor- 
mation available and the accuracy of sojsie of 
the analytical results. 

An irapingeirient -cooled blade was analyzed 
for an assumed cycle between effective gas tem- 
peratures at midspan of 1089 K (iSSgo F) to 
1839 K (2850° F) with a coolant temperature of 
533 K (500^ F), The cycle consisted of accel- 
eration and deceleration transients of 4 sec- 
onds each and a 30 -minute hold t-jjne at the high- 
temperatiu^e part of the cycle. Transient metal 
temperatures were computed from a quasi -three- 
dimensional heat transfer analysis program. 

The thermal transients and chordwise temperature 
gradients were subsequently altered to increase 
the severity of the thermal loading, and the 
mechanical loads were removed to simulate a vane 
problem. Inelastic strains were computed from 
a three -dijtiensional, finite -element computer 
program for the solution of nonlinear problems 
involving creep and plasticity (7) and (s)<, 
Cyclic life was calculated from Strainrange 
Partitioning relations for the reversed portion 
of the strain cycle and by time fractions for 
the unreversed steady-state cregr portion. In 
all, the design case and tliree parametric cases 
were analysed for the impingement -cooled blade 
configm’ation, 

Strainrange Partitioning was also applied 
to the inelastic strain liysteresis loops cam- 
puted for the design case of a film-impingement- 
Gooled vane in (4) by a one -dimensional stress 
analysis progi'sm. The predicted lives at crit- 
ical locations of this vane were ccaapared with 
experimental results from cyclic tests in a 
static cascade facility at the NASA Lewis Re- 
search Center reported in (4). 
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AmYTXGAL PROCEDUBE 


The heat transfer, stress analysis, and 
life prediction methods used for the 
iiiipingeiiient“Cooled blade study are described. 

The strain hysteresis loops for the filjn- 
ijnpingement-cooled vane were computed by using 
e:qperiinental etnd calculated transient temper- 
atures and a one -dimensional cyclic stress an- 
alysis program based on a kinematic hardening 
rule. These ai’e described in (4) and will not 
be further discussed herein. 

TRANSIENT TEMPERATURE ANALYSIS - A quasi- 
three -dimensional, time-dependent, thermal an- 
alysis computer program (9) was used to calcu- 
late the transient temperature distribution for 
the cyclic stress -strain analysis. The program 
computes the temperature and coolant flow dis- 
tribution for a turbine blade cooled by impinge- 
ment and crcssfloir convection. 

The blade selected for this analysis was 
a cast shell airfoil with an impingement insert. 
Eig. 1 shoT!7s a schematic cross-sectional view of 
the blade. Chord length and blade span are both 
3,8 cm (1.5 in.), and the wall thickness is tap- 
ered frojii 0.127 cm (0,050 in.) at the hub to 
0.076 cm (0.030 in.) at the tip. The impinge- 
ment insert is also tapered to maintain a uni- 
form coolant -channel width. The forvyard region 
of the blade is cooled by impingement, and the 
spent impingement air flows chordvfise through 
the trailing edge, cooling that region by forced 
convection. 

Input to the prograju includes a description 
of the geometry, an estimate of total coolant 
flow, the coolant supply pressure and temper- 
ature at the hub end of the insert, the hot-gas- 
side heat transfer coefficient distribution, the 
gas total temperatm’e distribution, and the gas 
static pressure at the trailing edge. For a 
transient calculation, tables of coolant supply 
pressure, hot-gas temperatures, hot-gas static 
pressuT’e at the trailing edge, and wheel speed 
as a function of time are required. 

The blade to be analyzed was divided into 
radial sections by chordwise planes and the lay- 
ers were treated individually, except that 
thermal communication between sections was main- 
tained in the three-dimensional heat conduction 
equations. The coolant -channel floiv for each 
layer was treated as one-dimensional compress- 
ible flow with friction, heat transfer, and mass 
addition. Fig. 1 shows how a layer is divided 
into calculation stations and the arrangement 
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of nodiSs at a station. 

The calculational procedure requires that 
on initial steady-state solution be determined. 
The calculations begin at the hub end of the 
blade and are marched radially through the 
blade, layer by layer, At each layer, the mean 
pressure and temperature in the coolant insert 
is calculated, assimiing adiabatic flov? with 
friction. These conditions are then used, along 
with the latest coolant -ehanne], pressure dis- 
tribution, to calculate the impingement jet 
flows and heat transfer coefficients. Using 
the new heat transfer data, the wall conduction 
finite difference equations and the coolant - 
channel energy equations are solved for a new 
temperature iistribution, and the coolant - 
channel momf;ntum equation is solved for a new 
pressure dlstributioni This is then used to 
recalculate Impingement flows, and the calcu- 
lations are repeated until foui* successive 
pressure distributions match within a specified 
tolerance band. This procedure is repeated for 
each layer of the blade. For the initial steady 
state, after all layers are coiaplete, the total 
coolant used is compared with the initial es- 
timate. If they do not natch within toleraxices, 
the inlet flow is adjusted and the calculations 
are repeated. Once the steady state has been 
established as an initial condition, the tran- 
sient calculations are started. The progrma 
solves the transient finite difference eqmtions 
for conduction in the blade wall and the time- 
dependent energy and momentum equations in the 
coolant channel. 

STKE03 AKALYSIS f-IETHOD - Analyses were 
performed ■5Jith a modified version of the Bias 55 
computer program for the in-core solution of 
structural problems involving time -independent 
inelasticity and thermally activated, time- 
dependent inelasticity. The theoretical back- 
ground and capabilities of this program are 
discussed in (7) and (8), The itiain plasticity 
and creep equations ore described in the appen- 
dix. Elas 55 was extended to accommodate prob- 
lems involving rotating machinery and reversed 
inelastic strain. 

Plasticity computations are based on the 
von Mises yield Condition and associated flox^ 
rule and assume isotropic hardening. The or- 
iginal version of Bias 55 used a viscoelastic 
Kelvin model to compute creep strains. However, 
this model was found to be inadequate to cor- 
relate the creep characteristics of the nickel- 
base superalloy blade material, cast IH 100. 
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To ittorease the accuracy of the creep computa- 
tioiis, a secondary creep rate term ms added to 
the creep eq,uations and the Kelvin model was 
used essentially to calculate the priinary creep 
strains. 

To avoid Introducing instahilitles, the 
creep ccmputations for steady-state hold times 
have to be performed in fine tiime increments. 
However, this increjaental approach reduces the 
accuracy of the calculated creep strains because 
the correlated creep characteristics are being 
extrapolated to strain lev.'^ls several magnitudes 
smller than is warranted b^' the experimental 
data. 

The analytical results were related to uni- 
ajcial stress-strain and creep data by combining 
the stress and strain components in terras of ef- 
fective (von Mises) stresses and strains. De- 
termining the sigris of the effective plastic and 
creep strains presented a problem since both 
eq.uations (A2) and (A12), in the appendix, al- 
ways give positive values. Therefore, signs 
were allocated to the effective strain incre- 
ments based on the signs of the normal stresses 
with the greatest magnitude during the tiitie in- 
crement. 

The program calcuiates average stresses and 
strains at the centroid of each element. Ther- 
jital load vectors were based on nodal point tem- 
peratures obtained frora the heat transfer anal- 
ysis. 

STjRESS -STRAIN DIAGl^Ml - The stress-strain 
curve for a virgin material during the first cy- 
cle is not representative of the stress strain 
behavior for subsequent cycles and may give rise 
to substantial errors if used for cyclic life 
prediction. Extensive cyclic tests of blade and 
vane materials have shown that, for a given 
strain range, the cyclic stress-strain hystere- 
sis loop will attain a stable form after a rel- 
atively small number of cycles, as typified by 
the loops shox-m in Fig. 2. 

Neither isotropic nor kinematic hardening 
specify this type of stress -strain behavior; in 
fact, neither of these hardening models will 
give a stable cyclic stress -strain diagram un- 
less the strain hardening parts of the diagram 
are straight and parallel lines. The type of 
stress -strain response exhibited in Fig, 2 
could be more accurately represented by i4x‘oz or 
mechanical sublayer models, as discussed in 
(lO). xlie Hiajor drawback in using either the 
toos or mechanical sublayer models in lar’ge- 
scale problems such as cyclic loading of tur- 
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bine "blades is fchs.t the storage requirements for 
these models are much greater than those for 
isotropic hardening. It is of interest that 
(10 ) also dejttonstrates that kinematic hardening 
can give less accurate results than isotropic 
hardening under multiajcial loading. 

To represent the actual cyclic stress - 
strain diagrajn by an isotropic hardening model, 
a rectangular diagratti with the same plastic 
strain range was assumed, as shcn-m in Pig, 2, 

The yield strength was selected so that the 
plastic strain energies represented by the en- 
closed areas of the isotropic and actual stress- 
strain diagrams were kept approximately the 
saiiie. The yield strengths are both temperature 
and strainronge dependent. Temperature depend- 
ence effects were obtained from cast IN 100 
monotonic stress-strain curves presented in 
(n). Strainrange effects were determined from 
the results of cyclic stress -strain tests at 
1200 K (1700® F) of cast IN 100 thin -wall tubu- 
lar speciifiens performed at the Lewis Research 
Center by using the techniques described in (12). 

These data are plotted in Fig. 3 in terras of 
stress and plastic strain ranges with zero mean 
stresses and strains. The cyclic frequency was 
0.5 Hz. 

To apply the strainrange corrections, the 
problem had to be run first with estlsnated 
yield stresses to obtain a preliminary solution. 

The yield strengths at the elements that ex- 
hibited srielding were then corrected on the ba- 
sis of Fig. 3 and the problem was rerun. Be- 
cause of the computing time involved, it was not 
possible to carry this process through to a com- 
pletely iterated solution. If the stresses at 
the most critical locations agreed reasonably 
well with the data of Fig. 3, the iteration was 
stopped at the second solution. In no case was 
the problem carried beyond a third iteration. 

STRESS ANALYSIS PROCEDURE - The airfoil 
shell was modeled by using 165 eight -node hexa- 
hedral elements in five spanwise layers with a 
total of about 1200 unsuppressed degrees of 
freedom; the modeling aroiind the airfoil cross- 
section is illustrated in Fig, 4. Blade taper 
and tilt were accounted for in deterifdning the 
cooi’dinates of the mesh network, 

The load-time history was set up so that a 
cycle consisted of 24 increments of temperatures Kaixfttian 
and mechanical loads; the latter included both 
centrifugal and gas pressure loading. A hold 
time of 30 minutes was applied at the high-gas- y 
temperatxire part of the cycle. To obtain a sol- 
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ution for a cycle, the problem ws run foi* an 
additional half cycle and the etrain cycle or 
hji-»teresis loop used for life prediction vjas 
determined from the inelastic strains after the 
first half-cycle. The program vjas too costly 
in computing time to extend the cycling beyond 
1 “ cycles. 

The program was stopped and restarted 20 
times for each case to change loads, temper- 
atures, time increments, and material proper- 
ties and to store the accumulated solution in 
case of a breakdown of the coBiputer system. 

About 200 000 words of core memory on the 
Univac 1110 computer were required to nm the 
problem. Each solution required approxlJitately 
5 to 7 hours of computing time. 

LISE PRBDICTIOH - Cyclic lives to crack in- 
itiation were calculated by losing Strainronge 
Partitioning when completely reversed inelastic 
strain cycles were obtained. For the unreversed 
strains (ratcheting), a linear-exhaustion-of- 
ductility concept should be used. In the pres- 
ent case, the only significant ratchet strains 
were creep strains induced diu'ing the 30 -minute 
steady-state hold time. Therefore, the amount 
of damage due to creep ductility exhaustion was 
based on time fractions (time of stress to rup- 
ture at saine stress) and the isionotonic creep 
rupeure data in (ll). This method was simpler 
and probably more accurate than using ductility 
exliaustion based on the creep strain increment 
during the first cycle. 

The Strainrange Partitioning relation pre- 
sented in Fig. 5 for cast IN 100 alloy were 
used for the isupingement -cooled blade. The cyc- 
lic stress-strain behavior shown in Fig. 6 and 
the Stralm'ange Partitioning relation in Fig. 7 
for cast Mar M 302 alloy were used in the com- 
putations for the filja-impingeiiient- cooled vane. 
These data were obtained from cyclic tests of 
tub'ilar specimens prepared froia the saJiie mater- 
ial billet as was used to fabricate the vanes, 
Bata fi'om tests conducted in a vacuum are re- 
ported in (13). Additional data obtained frowi 
tests conducted in air are also shown in Fig. 7,* 
these show considerably lower lives than the 
data of (13) and were, therefore, used for the 
vane life prediction, m (2) it is demonstrated 
that the Strainrange Partitioning relations are 
relatively insensitive to temperature for ma- 
terials ^jhose ductility is nou a strong fimction 
of temperature^ hovfever, such is not the case 
for Mar M 302. To make the Strainrange Par- 
titioning relations presented in Fig. 7 usable, 
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they v;ere corrected by the procedure described 
in (14), This procedure consists of multiplying 
the strain ranges by the ratio of the ductility 
at the temperature where inelastic straining in- 
itially occurred during the cycle to the ductil- 
ity at the uniaxial specimen test temperature, 

A probleia arises in applying the one- 
dimensional stress analysis results for the 
film-iittpingement- cooled vane to life prediction 
at stress risers since the analysis only com- 
putes nominal stresses and strains. The nominal 
strain ranges have to be multiplied by strain 
concentration factors at nodes ad^iacent to dis- 
continuities, such as the filjft-cooling slots and 
the inside wall in the leading-edge region where 
chordwise fins were present, as illustrated in 
Fig, 8, These strain concentration factors in 
the presence of inelastic strains 7?ere calculat- 
ed by the Neipber method described in (15) that 
used K- « k2/k . where K. is the strain con- 

centration factor, Kj. is the theoretical stress 
concentration factor and K is the ratio of 

tj 

the actual to noBiinal stress range. 

RESULTS Mm DISCUSSION - IMPINGEMENT -COOLED 
BLADE 


I-ffiTAL TEMPERATURE GRADIENTS - The chordwise 
temperature distributions and the transient 
thermal responses at the leading and trailing 
edges are summarized in Figs. 9 and 10, respec- 
tively, for 'C'he four conditions that were con- 
sidered: a design condition and three parajjiet- 

ric variations. 

Case 1 was the design condition. Temper- 
atures at the leading edge ^;ere relatively low, 
as shorn in Fig. 9(a); the hottest temperatures 
occurred in the trailing-edge region. The maj:- 
imujii chordwise teiaperature difference at mid- 
span was about HU K (200° F), During the ac- 
celeration and deceleration transients, the 
l.eading-edge temperature was lower than the av- 
erage midspan temperature (Fig. 10(a)), indicat- 
ing that the thermal stresses at the leading 
edge would ali.-ays be tensile. 

In case 2 the severity of the thermal load- 
ing was increased by increasing the steady-state 
temperature by an arbitrary ajnount of about 
Ul K (200° F) at the leading edge and decreas- 
ing it on the pressure side between the mid- 
chord and the,- leading edge, as shoim in 
Pig. 9(b). This doubled the maxijtiura chordwise 
temperattuce difference at midspan to about 222 K 
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(400° F). The transient response for case 2 vias 
also made more severe bjr reducing the leading- 
edge temperature about lH K (200° F) at the 
IcKf-temperatxire part of the cycle. As a result 
of these changes, the leading-edge temperature 
curve in Fig. 10(b) crosses the average temper- 
ature curve, resulting in compressive thermal 
stresses -,t the steady-state condj.tlon and ten- 
sile therjnal stresses at the beginning and end 
of the cj/'cle, 

A third condition (case 3) was analyzed 
with the same thennal load conditions as case 2 
but with fill the centrifugal and gas pressure 
loading removed. This changed the problem co 
one of a cantilevered vane. Since vanes are ex- 
posed to higher gas temperatures and, therefore, 
have to withstand higher metal temperatures than 
blades, a fourth condition (case 4) was consid- 
ered in which all the metal temperatures were 
raised 56 K (100° F) from case 3. Figures 9(c) 
and 10 (c) show the airfoil temperature distribu- 
tions and transients for ae.se 4. il.1 centrifugal 
and gas pressure loads were also removed for 
case 4 to simulate a vane, 

STRAIN CYCLES - The Computed strain cycj.es 
for cases 1 to 4 at three airfoil locations are 
presented in Fig, 11. Locations A, B, and G 
were selected as generally the most critical lo- 
cations in the leading-edge, tr ailing-edge, and 
pressure-side midchord regions, respectively. 

The most critical span position was usually at 
30 percent of the span length as measured from 
the airfoil base. In these cycles, almost all 
the inelastic strain at the steady-state condi- 
tion was due to creep, but plastic flov; predoasi- 
inated in the transient parts of the cycle. The 
computations shovred that the transient times 
were too short to induce significant creep 
strains. 

There was no inelastic strain reversal in 
case 1. Since the leading-edge thermal stresses 
in case 1 were tensile, and therefore additive 
to the centrifugal stresses, the leading-edge 
location A in Fig. 11(a) exhibited tensile creep. 

In all four cases, location B was a hot spot and 
location C was a cold spotj therefore, the creep 
strains were alvjays compressive at location B 
and always tensile at location C in cases 1 
to 4. Compressive creep strains were not fudged 
to be a danger to the structural integrity with- Kaufman 
out inelastic strain reversal and were not con- 
sidered in determining airfoil life, although 
they might eventually cause trailing-edge bm/ing 
in vanes. 
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When the l^?*Aing-edge tewperature was in- 
creased for cases 2 to 4, compressive creep oo- 
cm*red at location A. The only- condition imder 
which there was completely reversed inelastic 
straining was in case 2 where the compressive 
creep at location A was reversed by tensile 
plastic strain during deceleration (Fig. 11(b)), 
The loi’gest tensile creep strains also occurred 
in case 8. 

A comparison of Fig, 11(c) and (d) shows 
that the irialn effect of Increasing the metal 
temperatures 56 K (100^ F) was to cause a large 
increase in the plastic strain levels for the 
hot leading and trailing edges. Despite the 
higher temperature s, there t/as little change in 
the creep strain increments partly because of 
the lo^*7er yield stresses and, therefore, lower 
stress levels that resultedi Hotfever, the ac- 
curacies of the creep strain computations are 
problematic because of the extrapolation of the 
correlated creep characteristics to extremely 
small strains and time increments, as was dis- 
cussed in the STBESS AMLYSIS section. For this 
reason the relajced stresses at the end of the 
hold tixsje were considered to be a Kiore reliable 
basis on which to determine life than the ratch- 
eted creep strain increments. 

LIFE “ Ei’edicted cyclic lives for cases 1 
to 4 ore summarized in table l(a). All the 
computed lives were determined from ruptui'e due 
to tensile ireep by time fractions, except for 
the leading edge in case 8, where Strainrange 
Partitioning was applied. In all fooi’ eases the 
most critical location vjas at or near location 
G in Fig. 11, which approximately coincided 
with the lovjest temperature chordwise position 
in Fig, 9(b) to (d). 

In case 1, the predicted blade life was 
480 cycles at location G, The next most criti- 
cal location was at the leading edge, where the 
life would be 640 cycles. 

The only significant strain reversal oc- 
curred in case 2 at the leading edge, However, 
bhe predicted cyclic life based on the Strain- 
range Partitioning relation for the plastic 
tensile-compressive creep strain cycle is over 
100 000 cycles for the inelastic strain range 
at location A (Fig, 11(b)). Thus, increasing 
the leading-edge temperatur'e resulted in a hrige 
improveinent in leading-edge life for this con- 
fig iration. This is in line i-dth the temper- 
ature phasing benefits discussed in (14), which 
suggested that it might be beneficial to force 
the leading edge into compressive creep by 
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proper phasing of the temperatures at the edge 
and interior of the blade. The overall predict- 
ed improvement in blade life indicated in 
table 1(a) for ease 2 (760 cycles as compared 
with 480 cycles for case l) may be diie to the 
redu tion of about 17 K (30° F) in the average 
steady-state temperature shovm by comparing 
Figs, 10(a) and (b). Increasing the severity of 
the temperature gradient and transient response 
did not appear to deleteriously affect the blade 
life, these predictions ar-e in tj,ualitative 
agreement with results of cyclic tests of 
iidpingement -cooled airfoil specimens in the 
theriEial fatigue facility described in (16), 
these tests disclosed that a thin-wall, hollow 
airfoil with no stress concentrations, no webs 
to act as heat sinks, and no steady-state hold 
time was almost imjtiune to low-cycle fatigue 
failure even with leading-edge temperatures as 
high as 1311 K (1900° F) and chordwise temper- 
ature gradients as higli as 556 K (lOOO^ F). 

T'le predicted results for cases 3 and 4 are 
wh-' ';.f/,ld be intuitively e^cpected, Eemoving 
blade mechanical loads and reducing the 
problem to that of a cantilevered vane in case 3 
resulted in a substantial improvement in the 
airfoil life. Increasing the temperature levels 
by 5G K (100° F) in case 4 resulted in a reduc- 
tion in life front 2200 cycles for case 3 to 
1000 cycles for case 4. 

EESWfS Alffl DISCUSSIOH - F31M-IMPIKGEI4Effl- 
COOLED VAHE 


IffiTAL TEtffiEEATURE OMBIEi-JIS - Transient 
temperatiire distributions shown in Fig, IE were 
obtained from experimental data and theoretical 
heat transfer analyses for cycles over an effec- 
tive gas temperature range of 922 to 1644 K 
(1200° to 2500° f) with a coolant temperature of 
811 K (1000° F) , The vane hot spots were at the 
upstream sides of the pressure and suction sur- 
face film-cooling slots (locations B and 0), 

The leading-edge stagnation point was hotter 
than the bulk temperature of the vane during 
the acceleration part of the cycle and colder 
during deceleration. The inside wall at loca- 
tion A in the leading-edge region was the cold 
spot during most of the deceleration, 

STKAIH ClCiiES - The cohiputed noEtiirial , 
inelastic -strain versus temperature loops at the 
most critical locations are Ehcsm in Fig, 13 for 
the fifth cycle. The inelastic strains were 
predominately plastic since the transient times 
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were very short and the oyoles did not include 
any steady-state hold tiities. These loops were 
obtained from the nominal total strain cycles 
presented in (4) by subtracting the elastic 
strain components. 

The vane airfoil vfas modeled for the stress 
analysis with four layers of elements across 
the wall. Location A in the leading-edge re- 
gion was at the inmost element and adjacent to 
the chordwise fins. The fins could not be con- 
sidered in the one -dimensional stress analysis 
and, therefore, the computed inelastic strain 
range for location A was multiplied by a strain 
concentration factor for life prediction. Lo- 
cation B* in Fig. 13 was on the downstream side 
of the pres sure -surface, film-cooling slot and 
should not be confused with location B in 
Fig. IS, which is a hot spot on the upstream 
side. The nominal inelastic strain range for 
location B' in Fig, 13 did not tahe into account 
the discon' Inuity due to the adjacent film- 
coding slot and, therefore, was also multiplied 
by a strain concentration factor. By Keuber 
methods (15), the strain concentration factors 
were cdculated as 3,91 for the chordwise fins 
at location A and 5.59 for the corners of the 
filiifi -cooling slots at location B' . 

LIFE - The predicted lives presented in 
table 1(b) were obtained from the Strainrange 
Partitioning relation shown in Fig. 7 for ten- 
cile plasticity reversed by compressive plastic- 
ity for the vane alloy Mar M 302 in an air en- 
vironment. oince the ductility of this material 
is influenced by temperature in a kno^«m manner, 
the inelastic strain ranges presented in Fig. 7 
as a function of cyclic life were corrected, ac- 
cording to the effect of temperature on ductil- 
ity, by using the method of (14). A factor of 
0.477 for location A and a factor of 0.293 for 
location B' were applied to the strain ranges. 
These correction factors vjere based on the ratio 
of the ductilities at the temperatures where 
tensile inelastic straining initially occurred 
(1033 K (1400° F) at location A and 839 K 
(1050° F) at location B' ) to the ductility at 
the test temperature of 1273 K (1832° F) , 

The predicted life for complete fracture of 
60 cycles at location A is somewhat greater than 
the crack initiation life indicated by the ex- 
perimental evidence. At the end of 60 cycles 
(20 at the analytical condition and 40 at less- 
severe conditions) one of the test vanes was 
sectioned for inspection of the inside surface 
the leading edge. This inspection revealed 
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cracking at the base of many of the chordwise 
fins (location A), as shown by the photomicro- 
graphs of Fig, 14, The wide crack seen in 
Fig, 14(a) almost propagated to the outside sm— 
face; other signs of crack initiation ore vis- 
ible in Fig, 14(b), The greater predicted life 
for location A may be due to the fact that the 
Mar M 302 test specimens used to obtain the data 
in Fig, 7 were cast smooth, machined, and pol- 
ished; the internal surfaces of the vane wo\£Ld 
have much rougher surfaces and, therefore, low- 
er lives. 

The predicted life of eight cycles at loca- 
tion B' is consistent with the e:sperimental re- 
sults, although the exact ^lumber of cycles at 
which failure began is uuL'fnovm. Fatigue cracks 
were found on the downstreain side of the 
pre s sure -svir face, film-cooling slots after the 
60 cycles of testing. Similar cracks were found 
in 69 out of 72 vanes of this configuration 
(Fig, 15), many of which had undergone less than 
40 cycles in a research turbojet engine operat- 
ing at essentially the some conditions. The 
leading-edge crack shw-m in Fig. 15 was revealed 
by metallographic inspection to have initiated 
at the inside surface and propagated through 
the wall. 

The severity of the fatigue problem for 
the film-impingement -cooled airfoil configura- 
tion was caused by the presence of the web, 
which remained cold during acceleration and hot 
during deceleration and constrained the move- 
ment of the airfoil walls. Even if the effects 
of the stress risers were neglected, uhe in- 
elastic strain ranges for locations B' and A 
were large enough that failures would have been 
predicted in 53 and 270 cycles, respectively. 

The thermal gradients and levels shown for the 
film-impingement -cooled vane in Fig. 12 appear 
to be no more severe than those shown in Figs. 
9(c) and 10 (c) for case 4 of the impingement - 
cooled airfoil. 

SUt#IARY OF RESULTS 

The inelastic strain cycles and the cyclic 
lives to crack initiation of an Impingement - 
cooled blade under various thermal and mechanic- 
al loads and of a filiri-impingement -cooled vane 
were calculated from cyclic stress analysis pro- 
grams, using computed transient metal temper- 
a.5\ires and either a Strainrange Partitioning or 
time fractions approach, where each was applic- 
able. The results of these analytical studies 
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can be simimarized as follows: 

1. The life predictions made by using 
Strainrattge Partitioning and time fraction meth- 
ods were generally consistent with, although 
not completely verified by, limited experimental 
evidence from cyclic tests in research cascade 
facilities of the saiiie or similar cooling con- 
figurations . 

2. The Strainrange Partitioning approach 
appears to be applicable to a constrained air- 
foil construction, such as in the filiti- 
impingement -cooled vane. In a constrained air- 
foil it is relatively easy to induce lax'ge in- 
elastic strains and strain reversal during the 
transient parts of the cycle. 

3. Time fractions were more applicable for 
life prediction with the f'ully impingement - 
cooled blade than Strainrange Partitioning. In 
a thin-wall, hcllow airfoil construction with 
no webs, such as exists in the ijiapingement- 
cooled blade, it is difficult to induce inelas- 
tic strain reversal even \inder severe thermal 
transients. The dmaage due to inelastic strain 
hysteresis was of minor consequence compared 
with the monotonic creep ratcheting damage. 

4. Increasing the severity of the thermal 
transients and gradients for the impingement- 
cooled blade did not deleteriously affect the 
predicted lives.- Raising the leading-edge tem- 
perature above the average span-section temper- 
ature actually increased the predicted life of 
the leading edge by putting it in compression. 

5. As expected, removing the airfoil mech- 
anicii?. loads to siinulate a vane while maintain- 
ing the same thermal loads increased the pre- 
dicted life. However, changing the thermal 
loads by increasing the metal teraperatui’e levels 
56 K (100° F) resulted in a reduction in pre- 
dicted airfoil life of about 50 percent. 

APPENDIX - ELASTICITY-CREEP STRESS AND STRAIN 
EQUATIONS 


Effective stress is represented by 


a “ 



S S. . 
10 


(Al) 


and the effective plastic strain increment is 
represented by 


d€^ — 
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where incremental plastic strain 

vector and is the stress-flevlation tensor. 

The components of are 


S = CT “ 
X 'X 


fo + a + 0 ) 


X 


(AS) 


S 


y 


y 


s 


z 


( + gy t 

3 

( g^ + Oy + g^) 

3 


(A4) 

(AS) 


vjhere o' ,g , and g denote the noritial stress 
3c’ y z 

components in the directions of the global axes 
and 
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T 

denoting the shear 



xy xz zy 

components. 

The uniaxial secondary creep rate is ex- 
pressed in exponential form 


ic 


= Ag^ 


(A9) 


where A is a constant and n is a fvinction of 
temperatiire. 

The multiaxial secondary creep rate is 


G . ~ £ A g S . 


(AlO) 

id 2 

At the end of time t the creep strain is 


= ^i.i 

E., 

ij 


U fi - 7 (All) 

J J L 


where a and e ar-e the stress and creep 
ij ij 
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axe the 


strain vectors and and ti?. 

ij Xj 

Ifelvin model spring and dashpot constants, re- 
spectively. 

The effective incremental creep strain is 
given by 


= y| 


(A12) 


which is analogous to equation (A2) for the in- 
cremental plastic strain. 
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Table 1. - Predicted Cyclic Lives 


(a) Iiftplngemejit -cooled blade 

Leading-edge Pres sure -surface Other critical 
Case location A location G loca t ions 


Number of cycles to failure 


1 640 

2 ^>100 000 

3 

4 


480 

1700 

2200 

1000 


ayeo 


(b) Film-impingement -cooled vane 

Leading-edge Pressure -side 
location A location B' 


Number of cycles to failure 


-60 


^Pressure surface ad;jacent to location C. Loca- 
tions A and C are at 30 percent of span and are 
shown in Fig. 11, 

“Strainrange Partitioning used. Other lives were 
based on time fraction, due to creep. 

*^Strainrange Partitioning used. Locations A oad 
B are sho-^-m. in Fig, 13. 
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Fiqure 1. * Cross section of impingement -cooled blade. 


STRESS -STRAIN DIAGRAM 
ISOTROPIC HARDENING MODEL 


STRESS 

t 



CURVE 1 


PLASTIC 

strain 




300 

- 200 

a 






< 

“ X 

»/» 

•^3 


2000 

1000 
wo 
600 

voS^ 

10 ' 


Oo • I MtelO^ UUpA^ INSI UNITS; 
Ao •0./00»l0‘u^pA*^IN 
CUSTOMARY UNITS 



J Lj. i l 111 J J l_J_ij.l-J I 


10* 

PLASTIC STRAIN RANCl. U. 


10 ‘? 2 , 10 ? 


cm/cm 


Tlqurt ). - Cyclic tlrm-slrtin twhavior ol cJtt IN 100 it 1200 K iirao^Fi 
Friqutncy, 0. S Ki; &\Np • 0. S7K Ao Ocp 



Figure 4. - Finite-element modeling of airfoil cross section for 
impingement -cooled blade. Airfoil chord length, 3. 81 cm U. SO in. I; 
airfoil span length, 3. 81 cm II. 50 in. ). 
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Figure 7. - Partitioned strain range versus life relations 
for Mar M 30? alloy at 1273 K 11832° FI. 
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Figure 8, - Film-impingement -cooled vene. 
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Figure 9. - Chordwise temperature distribution at high -gas- 
temperature, steady-state condition for impingement- 
cooled blade. 
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Figure 10. - Metal temperature cycle at midspan (or 
impingement -cooled blade 
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Fiqure 11. - Strain -temperature cycles for impingement -cooled blade 
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Figure 12. * Metal temperature cycle at midspan for film* 
impingement-cooled vane. 
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Figure 13. - Strain -temperature loops for film 
impingement-cooled vane. 




(a) CRACK PROPAGATION THROUGH AIRFOIL WALL 




(b) CRACK INITIATION AT BASE OF FINS 

Figure 14. • Fatigue cracks at inside surface of leading edge of film-impingement- 
cooled vane 
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Figure 15. - Fatigue cracks in film-impingement-cooled vane. 


